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Determination of the precise location of the multicritical point and phase boundary is a target 
of active current research in the theory of spin glasses. ^"^^-^ In this short note we develop a duality 
argument to predict the location of the multicritical point and the shape of the phase boundary in 
models of spin glasses on the square lattice. 

The first system we treat is a random Zq model with gauge symmetry which includes the ± J 
Ising model and the Potts gauge glass. Following the notation of Wu and Wang,^-^) the partition 
function is 

Z = 5]exp|^y(Ci-0 + Ji,)i. (1) 

Here, ^i{= 0, 1, • • • , g — 1) is the g-state spin variable, Jij{= 0, 1, • • • , g — 1) is the quenched random- 
ness, and V{-) is a periodic function with period q. The sum in the exponent runs over neighbouring 
sites on the square lattice. It is straightforward to generalize the duality transformation in the Wu- 
Wang formalism^^^ to the random model, and the result is 

Z = ^'exp| j;y(r,,, + Ji,)i (2) 
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where rjij denotes Xij is the bond variable on the dual lattice corresponding to rjij and e 

is the Fourier transform of e^. The prime in the sum indicates that the bond variables rjij and Xij 
are constrained such that their sums over a closed loop vanish. Thus the duality transformation 
consists in replacing the exponential in (2) with that in (3): 
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Average over quenched randomness is dealt with by the rephca method. A generahzation of 
the existing argument-*^^) suggests that it is convenient to consider the following quantity 

q-l (q-1 

" = '—' (5) 
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where pi is the probability that Jij takes the value I. It is not difficult to check that the duality 
transformation (4) changes a into a = (/"/a. The subspace defined by a = q^/'^ in the parameter 
space {pijP2j - ■ ■ ,T,- ■ ■) is therefore left invariant by the duality transformation: A point in this 
subspace is transformed into another point in the same subspace although the latter point does 
not necessarily have real values of Boltzmann weights. In general it is possible to consider complex 
values of Boltzmann weights (or parameters such as the temperature); the subspace a = g"/^ 
embedded in the complex- valued parameter space maps onto itself by the duality transformation (4). 
We shall consider the cross section of this subspace with the phase diagram with real parameters, 
which we shall call the invariant subspace in the following. 

One should be careful in identifying the invariant subspace with a critical surface. We discuss 
this problem later; let us see here what the invariant subspace looks like in the quenched system, 
which is obtained by taking the limit n ^ of the relation a = g"/^: 

E^^zlogfxJe^C'+O-no] =liog,. (6) 
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We now apply (6) to a few examples and discuss the consequences. The first case is the ±J 
Ising model {q = 2) with ^(0) = /3J, V(l) = —[3 J and po = p,pi = 1 — p. The invariant subspace 
(6) then reads 

p log(l + e-2/3-^) + (1 - p) log(l + e^^') = ^ log 2. (7) 

This curve is plotted in Fig. 1 in the p-T phase diagram. It is observed (and can be confirmed 
analytically) that the curve reaches the point with minimum p on the Nishimori line^^^ e~^^'^ = 
(1 — p)/p, from which and (7), the minimum point Pc is found to satisfy 

log 2 

-Pc \ogpc - (1 - Pc) log(l - Pc) = (8) 

or Pc = 0.889972. This coincides with numerical estimates of p at the multicritical point with high 
precision: 0.8905(5),^) 0.886(3),2) 0.8872(8),^) 0.8906(2)^) and 0.8907(2).^) 

The next example is the g-state Potts gauge glass with V{0) = PJ,V{1) = V{2) = ■■■ = 
V{q — 1) = and Po = ^ — {q — ^)p,Pi = ■ ■ ■ = Pq-i = p, where we have followed the notation of 
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Jacobsen and Picco.^^ The invariant subspace (6) is, in this case, 

{l-((7-lMlog{l + (g-l)e-^-^} 

+(g-l)plog(g-l + e/3'^) = i^. (9) 

This formula reduces to the Ising counterpart (7) when g = 2 under appropriate changes of energy 
scale and notation (exchange of p and 1 — p). The curve specified by (9) reaches the extremum 
point in the p-T phase diagram again on the Nishimori line p = l/(e^'^ + q — 1)-^^^ In particular, 
for q = 3, this extremum point Pc satisfies 

-(1 - 2p,) log(l - 2p,) - 2pe logpc = (10) 

whose numerical value Pc = 0.0797308 coincides very well with the location of the multicritical 
point, 0.079-0.080.6) 

It is possible to apply the same argument to the random Ising model with a general distribution 
function P{Jij) of exchange interactions.^^) The result for the invariant subspace is 

dJ,, P{J,j) log(l + e-'f^'^^) = (11) 

-oo ^ 

This formula applied to the Gaussian model with P{Jij) oc Q-iJa-Jo)'^ I'^J^ represents a curve in the 
phase diagram similar to that of the it J model in Fig. 1 and again has the smallest Jq/J on the 
Nishimori line /3J^ = Jq}^^ This point is at Jo/ J = 1.02177 which is very close to the multicritical 
point evaluated numerically (Horic, Nemoto, Hukushima and Ozeki, private communication). 

What do all these results mean? One possibility is that the exact locations of multicritical 
points have been derived for the above models. The present argument using the subspace invariant 
under duality is, in general, not guaranteed to give the exact phase boundary or multicritical point. 
Nevertheless we have several reasons to conjecture that our result may be exact, in particular for 
the multicritical point. 

The first reason is that the invariant subspace a = (pl'^ coincides with the exact phase boundary 
for the zb J Ising model in the case of n = 1 and n = 2.^^) It should be noticed that this coincidence 
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Fig. 1. Invariant subspace of the ±J Ising model is the curve drawn by the solid line. Shown dashed is the Nishimori 
line. 
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exists only above the multicritical point in the phase diagram for the n = 2 case, which may be 
related to the strange shape of the curve below the multicritical point in Fig. 1. The second reason 
is the impressive agreement with numerical estimates explained above, which seems to us beyond 
a simple coincidence. The third evidence is the fact that the Nishimori line appears naturally from 
the present duality formalism. It is well established that the multicritical point is located on the 
Nishimori line for models with gauge symmetry. The natural emergence of the Nishimori line 
from arguments without explicit use of gauge transformation strongly suggests that something deep 
may be hidden behind the scene. 

There are of course several problems that deserve special caution. As pointed out by Aharony 
and Stephen, duality does not yield fixed points of the transformation for random models whereas 
fixed points of duality are often identical to critical points in non-random systems. Aharony and 
Stephen thus argued that duality in random systems, unable to identify fixed points, is not a useful 
tool of analysis. We are suggesting here that an invariant subspace, if not fixed points, may be of 
some use in random systems. 

Another problem is whether or not the whole invariant subspace (6) or (11) coincides with the 
phase boundary. It is dangerous to accept this identification at least below the multicritical point 
in the phase diagram because of the itJ Ising model with n = 2 as mentioned above. The origin 
of disagreement of our result (7) near p = 1 with the perturbative calculation of Domany^^^ should 
also be clarified in future investigations. 
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